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We improved the decoupling approximation of the double-time Green's function theory, and applied 
it to study the spin-i two-dimensional antiferromagnetic Heisenberg model with broken bonds at 
finite temperature. Our decoupling approximation is applicable to the spin systems with spatial 
inhomogeneity, introduced by the local defects, over the whole temperature region. At low temper- 
atures, we observed that the quantum fluctuation is reduced in the neighborhood of broken bond, 
which is in agreement with previous theoretical expectations. At high temperatures our results 
showed that the quantum fluctuation close to the broken bond is enhanced. For the two parallel 
broken bonds cases, we found that there exists a repulsive interaction between the two parallel 
broken bonds at low temperatures. 



PACS Numbers: 75.fO.Nr, 75.30.Hx, 75.50.Cx 
I. INTRODUCTION 

In recent years, the two-dimensional (2D) antiferro- 
magnetic (AF) spin system on a square lattice has been 
one of the subjects of major interest in condensed matter 
physics 0]. This follows from the experiments that cop- 
per oxide sheets in the high-Tc superconductors show 
strong AF correlations [1). The undoped copper oxide 
materials are layered AF insulators and well described by 
the 2D AF Heisenberg model |Q] . Doping holes into these 
materials leads to frustration of spins and ultimately to 
destruction of AF long-range order (LRO). In the ex- 
treme limit of static hole, holes act as local defects and 
the inhomogeneous Heisenberg model is believed to de- 
scribe some of the physics ||^ . 

Several numerical and analytical works have been de- 
voted to the effects of isolate defects, e.g., static vacancies 
[Q-^ , broken or ferromagnetic bonds J/j-p^ , and dynamic 
holes , on the magnetic properties of the 2D antifer- 
romagnet. The inhomogeneous Heisenberg systems are 
mainly divided into two types. One is the site-defect 
(SD) model and the other one is the bond-defect (BD) 
model. These models are important for many fundamen- 
tal problems, such as frustration, phase separation, and 
spin glass. Here we adopt the BD model to study the 
effects of broken bonds replacing AF links in the spin-^ 
Heisenberg model at finite temperature. For zero tem- 
perature, Lee and Schlottmann have studied this model 
using the linear spin wave (LSW) theory §. They ob- 
served that the quantum fluctuation is reduced in the 
neighborhood of the impurity link and the local magnetic 
moment is enhanced, in agreement with results obtained 
by Bulut et al for static vacancy case |^ . 

At any finite temperature, Mermin and Wagner | [T^ ] 
have proved that for models such as the AF Heisenberg 
model the AF LRO is destroyed by strong thermal fluc- 
tuations in low dimensional systems. Spin wave theory 



which is based on the existence of LRO can not be di- 
rectly used at finite temperature. Alternatively, Kondo 
and Yamaji jl^ proposed a second-order Green's func- 
tion (SOGF) theory to study the low dimensional Heisen- 
berg model over the whole temperature region. At high 
temperatures, this theory reproduces the correct results 
obtained by the high temperature expansion method [ p^ . 
On the other hand, the results at low temperatures are 
similar to those of the modified spin- wave theory [ p^ . 
The SOGF theory does not violate the sum rule and 
rotation symmetry of spin correlations. In the SOGF 
theory, the decoupling approximation is at a stage one 
step further than Tyablikov's random-phase approxima- 
tion (RPA) The SOGF theory has been successfully 
applied to various low dimensional homogeneous systems 
without LRO, such as the one-dimensional (ID) Heisen- 
berg model pl^, the ID XXZ model the 2D Heisen- 
berg model and the 2D XXZ model Their re- 
sults are in qualitative agreement with those numerical 
results |20y2l|] over the whole temperature region. 

In this paper we extend the SOGF theory to study the 
inhomogeneous Heisenberg model by improving the de- 
coupling approximation. The decoupling approximation 
proposed by Kondo and Yamaji (KYDA) |13| can not 
be directly applied to the inhomogeneous case without 
modification. In order not to violate the sum rule in the 
inhomogeneous case, we introduce an improved decou- 
pling approximation, which is equivalent to the KYDA 
in the homogeneous case. In our approximation two pa- 
rameters Pi and (3j are attached to the correlation func- 
tion of the two corresponding spins on sites i and j as 
(3i{Sf S~)(3j. While in the KYDA, only one parameter 
ai is introduced. It is clear that these two parameters 
account for vertex correction of spin-spin correlation and 
have to be introduced in order not to violate the sum 
rule of correlation functions. Thus, in our decoupling 
scheme, iV vertex correction parameters {j3i) for a lattice 
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of A'' sites were introduced and no other extra parameter. 
In the homogeneous case, we obtain = a and our de- 
coupHng approximation reduces to the KYDA. Applying 
this method to the spin-^ 2D antifcrromagnetic Heisen- 
berg model with broken bonds, we obtain reasonable re- 
sults over the whole temperature region. 

The present paper is organized as follow: in Sec. II we 
present our extension of the SOGF theory to the inho- 
mogeneous spin systems and discuss the improved decou- 
pling scheme. Our numerical results for some particular 
configurations of one, two, and three broken bond cases 
are studied in Sec. III. Finally, we conclude our findings 
in Sec. IV. 

II. EQUATION OF MOTION AND DECOUPLING 
APPROXIMATION 

The 2D spin-i Heisenberg model with bond-dependent 
exchange constants can be expressed by the Hamiltonian 

H = Y. M\{StS- + S-S+) + S^S]}, (1) 

where (ij) denotes a sum over nearest- neighbor (NN) 
bonds, and Jy is the exchange interaction between spins 
on sites i and j. For the homogeneous case, Jij is equal 
to J for all bonds. For the inhomogeneous case we are 
studying here, Jij equals to zero for broken bond, while 
it equals to J for unbroken bonds. 

We define spin Green's functions G by 

Gi^~J, t~t')^{{S!it); 5|(0)>. (2) 

After time-Fourier transformation, the equation of mo- 
tion of the spin Green's function G can be evaluated as 

LoG{i-j, ^) = ^jM+,((s+5r+,-^i,^r; s])), (3) 

with r]=x, y. The SOGF appears on the right hand side 
of Eq.(3). Furthermore, establishing equation of motion 
of the SOGF, we have 

= 2{S^Si_^_,^){di+rij - Sij) + 2Jj^i+.^((5f - S-^jj-, Sj)) 

(4) 

and the third-order Green's functions appear on the 
right-hand side of Eq. (4). 



In the homogeneous case, Kondo and Yamaji decou- 
pled the chain of equation (4) in an approximate way 
|l3t , for example 

{{S+S-S!; a{S+S-){{S^; (5) 

for i ^ m ^ n. The parameter a was introduced in order 
not to violate the sum rule of correlation functions. For 
the inhomogeneous case, the lattice translational invari- 
ance does not exist and one needs to introduce N such pa- 
rameters for a lattice of N sites. However, simply replac- 
ing a by ai leads to difficulty in solving self-consistent 
equations of G. Instead, we introduce N parameters f3m 
for each site m according to the following relations: 

~ S^[/3jn — 2(1 — (3rn)S^] 

S~ = - 2(1 - P,n)S^]S- (6) 

for the spin= i case. The decoupling approximation for 
the inhomogeneous systems is thus expressed as 

{{S^Sj^S^; Sj)) — {{S^^[f3„i — 2{1 — (3,yi)S^] 

[/5„-2(l-/3„)5^]5-; S])) 
^/3„(5+5-)/3„((5f; (7) 

Here we only keep terms of the lowest order (three op- 
erators) Green's function. On the right hand side of Eq. 
(7) two parameters f3m and /3„ are attached to the corre- 
lation function of the two corresponding spins on sites m 
and n. It is clear that these parameters are vertex cor- 
rections of the spin-spin correlations. According to the 
definition, we introduce one vertex parameter for each 
site, and no other extra parameters. In the uniform case 
we obtain f]'^ — a, and our decoupling approximation is 
the same as the KYDA (as shown in Eq. (5)). 

With the help of the above decoupling scheme, we ob- 
tain the equations of motion of the spin Green's function 
G 

uj^G{i - j, = ^ 2 Ji^j++^{C(i, i + rj){5i+n. i ~ ^% j) 
+ ((2J.,.+,(5f -5f+,) + H; (8) 

where 

H = ^ {2J,,,+^/[/3,+^./3,+^C(i + ri, i + r]')S^ 

-AA+„C(^ + ^?,^)5f+^,] 
+2J. +ri—ri',i+ri [A+^/3,C(i -I- 77, i)S-^^_^, 

AC(* + r/, z)5f+,)]} , (9) 

here the relation G{i,j) = (S+S^) = 2{S^S]) has been 
used. 

Since there exists no translational invariance, we can 
only solve this set of equations of Green's function in 
real space. For a finite lattice, the Green's function G 
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can be expressed in a matrix form G, and Eq. (8) can be 
rewritten as 



hG = C, 



(10) 



where matrices h and C are made of the nearest- 
neighbor, the second nearest-neighbor, and the third 
nearest- neighbor spin-spin correlation functions. We 
solve self-consistent equation (10) to determine the spin- 
spin correlation functions and the vertex correction pa- 
rameters /3i for each site so to calculate thermodynamical 
quantities at any temperature. 
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and found that its leading finite-size dependence is of 
order 1/-L, which agrees with the scaling law proposed 
by other theories |^,^. The extrapolated estimate of 
the staggered magnetization for an infinite lattice is m = 
limi^oo \/S{Q) — 0.23, which grossly agrees with the 
results of the other theories |^,^ . Even though we are 
limited to finite size clusters, our calculation can give 
reasonable estimates for the infinite system for problems 
we are interested in. 

Let us study the one broken bond case. All con- 
figurations of one broken bonds are equivalent within 
the periodic boundary conditions. Our numerical re- 
sults of the NN correlation functions near the broken 
bond for the 8x8 lattice at temperature T ~ 0.01(2J) 
are shown in Fig. 1. Our results show that the bro- 
ken bond enhance correlation between spins close to it, 
which means that the quantum fluctuation close to the 
broken bond is reduced. The biggest NN correlation func- 
tion is C(2, 5) ~ —0.224, which is about 8% lower than 
Ci — —0.207 of the uniform case. At zero temperature, 
Lee and Schlottmann have studied this model by using 
the LSW theory. They showed that the quantum fluctu- 
ation is reduced in the neighborhood of the impurity link 
and the local magnetic moment is enhanced. Our results 
agree with the results of the LSW theory. 

The scaling effect of the NN correlation functions close 
to the broken bond is shown in Table L We found that 
the results of the 8x8 lattice are very close to that of 
the 10x10 lattice. We also found that the NN correla- 
tion functions in the next neighborhood of the broken 
bond are reduced, that is, the quantum fluctuation gets 
enhanced as the distance to the broken bond increases. 



FIG. 1. The NN correlation functions of the 8x8 lattice 
with one broken bond at temperature T = 0.01(2J). The 
dotted line between site 1 and 2 (open circles) represents the 
broken bond. 

III. NUMERICAL RESULTS 

We have performed numerical calculation on the 6x6, 
8x8, and 10x10 lattices with periodic boundary con- 
ditions. For the convenience of comparison, we first 
study the 6x6, 8x8, and 10x10 homogeneous Heisen- 
berg antiferromagnet. Our results of these frnite lattices 
are quite close to the results of an infinite lattice p8| , 
since the interaction is short-ranged. For the ground 
state we obtained that the NN correlation function Ci 
is equal to -0.2080 for 6x6 lattice, which is quite close to 
the corresponding result —0.2067 for the infinite lattice 
p8[ . We also studied the staggered structure factor S{Q) 
{Q = (7r,7r) is the AF momentum), defined by, 

^w) = ]^E([(-i)'^^^''^'(*-^^)]') ' (11) 



TABLE I. The scaling effect of the NN correlation functions 
of some special NN sites (labeled as in Fig. 1). The NN 
correlation function of uniform case is Ci = —0.207. 



TV x iV 


C(l,2) 


C(3,4) 


C(2,5) 


C(5,6) 


6x 6 
8x 8 
lOx 10 


-0.114 
-0.105 
-0.101 


-0.200 
-0.202 
-0.203 


-0.226 
-0.224 
-0.224 


-0.204 
-0.205 
-0.205 



The energy cost for removing one bond from the ho- 
mogeneous lattice is defined as 

Ai = El - El (12) 

where -Eq is the ground state energy of the uniform case, 
and E\ is the ground state energy of the one broken 
bond case. For the 8x8 lattice, we obtained that Ai = 
0.10(2 J), which is smaller than the energy per bond = 
— 0.312(2J) of the 8x8 homogeneous Heisenberg model. 
This different is due to the fact that the magnetic system 
has lowered its energy considerably by readjustment of its 
spin correlations. 
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function Ci of the homogeneous Heisenberg model is also 
plotted (dotted line). In the low temperature region, the 
correlation function C(2, 5) is larger than Ci. When the 
temperature increases, C(2,5) drops very quickly and 
becomes smaller than Ci as T > 1.25(2J). For the 
8x8 lattice, the NN correlation function near the broken 
bond at temperature T = 1.5(2J) are shown in Fig. 3. 
The biggest NN correlation function is C(5, 6) = —0.168, 
which is only about 1% lower than the NN correlation 
function of the homogeneous case (Ci = —0.166). It is 
obvious that the effect of broken bond on the nearby NN 
correlation functions in the high temperature region is 
weeker than that in the low temperature region. 



"'^0.0 0.5 1.0 1.5 2.0 

T(2J) 

FIG. 2. The temperature dependent of the NN correlation 
functions C(2, 5) (solid line), C(5, 6) (dashed line) and C(3, 4) 
(dash-dotted line) of 8x8 lattice. The dotted line is the result 
of the homogeneous model. 

The temperature dependent of C(2,5), C(5,6) and 
C(3, 4) are shown in Fig. 2. For the convenience of com- 
parison, the temperature dependent of NN correlation 




FIG. 3. The NN correlation functions of 8x8 lattice with 
one broken bond at temperature T = 1.5(2J). The dotted 
line between site 1 and 2 (open circles) represents the broken 
bond. 



FIG. 4. Two different configurations of the two broken 
bonds case. 4(a) the two broken bonds are parallel and adja- 
cent; 4(b) the two broken bonds are still parallel but have a 
distance of two lattice spacing. 

For the two broken bonds case, we mainly consider two 
configurations as shown in Fig. 4(a) and 4(b). In Fig. 
4(a), the two broken bonds are parallel and adjacent. 
While in Fig. 4(b) the two broken bonds are still parallel 
but the distance between these two bonds increases to 
two lattice spacing. We obtained that the energy cost 
of constructing the configuration 4(a) from the homoge- 
neous lattice, e.g., removing two adjacent AF bonds, is 
^2a = 0.31 (2 J), which is higher than the energy of in- 
troducing two isolated broken bonds 2Ai — 0.20(2 J). In 
addition, is quite smaller than the two bond energy 
2^]^ = — 0.624(2J) of the uniform case. The correspond- 
ing energy of the configuration 4(b) is A2f, = 0.24(2 J), 
which is more closer to 2Ai than A2a. Thus configura- 
tion 4(b) has energy lower than that of 4(a). Our results 
can be interpreted as that there exists effective repulsive 
interactions between the two parallel broken bonds. 
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FIG. 5. The spin staggered structure factor S{Q) as a 

function of temperature T for the homogeneous lattice (cir- 
cle), the one broken bond case (square), the two parallel bro- 
ken bonds as Fig. 4(a) (diamond) and Fig. 4(b) (up triangle), 
and the three parallel adjacent broken bonds case (inverse tri- 
angle) . 

The spin staggered structure factor S{Q) {Q = (tTjTt)) 
as functions of temperature T for the homogeneous lat- 
tice and some particular configurations of one, two and 
three broken bonds cases are plotted in Fig. 5, respec- 
tively. The configuration of three broken bonds shown in 
Fig. 5 is that the three broken bonds arc parallel and ad- 
jacent. Our results showed that, although the AF corre- 
lation functions close to the broken bonds are enhanced, 
the AF order of the whole system will be suppressed as 
number of broken bonds increase. 



IV. SUMMARY 

In conclusion, we have extended the second-order 
Green's function theory to the inhomogeneous Heisen- 
berg model by improving the decoupling approximation 
introduced by Kondo and Yamaji. The Kondo and Ya- 
maji's decoupling approximation can not be directly ex- 
panded to the inhomogeneous case. In this paper, we in- 
troduced an improved decoupling approximation, which 
is in accordance with the Kondo and Yamaji's decou- 
pling approximation in the homogeneous case. In our 
approximaton two parameters arc attached to the corre- 
lation function of the two corresponding spins, which are 
vertex correction of spin-spin correlation and they have 



to be introduced in order not to violate the sum rule 
of correlation functions. We have tried the one parame- 
ter scheme (Eq. (5)) and met difficulties in getting the 
self-consistent equation (10) converge. In our decoupling 
approximation, there are N vertex correction parameters 
for a lattice of A'' sites and no other extra parameter was 
introduced so we did not add one more parameter than 
the one parameter scheme. Moreover, our decoupling ap- 
proximation reduces to the Kondo and Yamaji's in the 
homogeneous case. 

We apply this theory to study the 2D Heisenberg an- 
tiferromagnet with broken bonds for all temperatures. 
At low temperatures, our numerical results showed that 
the AF nearest-neighbor correlation functions close to 
the broken bonds are enhanced. That is, the quantum 
fluctuation is reduced in the neighborhood of the broken 
bond at low temperatures. Our results are in agreement 
with the results of other theories. As the distance to the 
broken bond increases, the NN correlation functions de- 
creases. By contrast, at high temperatures our results 
showed that the quantum fluctuation close to the broken 
bond is enhanced. For the two broken bonds case, we 
found that there exists repulsive interaction between the 
two parallel neighbor bonds. 
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